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A recent paper by I Webh et al. (2014c) on multi-symplectic magnetohydrodynamics 


(MHD) using Clebsch variables in an Eulerian action principle with constraints is further 
extended. We relate a class of symplecticity conservation laws to a vorticity conservation 
law, and provide a corrected form of the Poincare-Cartan differenti al form formu l ation o f 
the system. We also correct some typographical errors (omissions) in Webb et al. I ( 2014cl) . 
We show that the vorticity-symplecticity conservation law, that arises as a compatibility 
condition on the system, expressed in terms of the Clebsch variables is equivalent to 
taking the curl of the conservation form of the MHD momentum equation. We use the 
Cartan-Poincare form to obtain a class of differential forms that represent the system 
using Cartan’s geometric theory of partial differential equations. 


1. Introduction 

Multi-symplectic equations for Hamiltonian systems with two or more independent 
variables have been developed as a useful extension of Hamiltonian systems with 
one evolution variable t. This development has con n ection s with dual variational for¬ 
mulations of traveling wave problems (e.g. Bridget 1 1992 1). and is useful in numeri¬ 
cal schemes for Hamiltonian systems. Bridges and co-workers used the multi-symplectic 
approach to study linear and nonlinear wave propagation, g eneraliza t ions of wave ac¬ 
tion, w av e modulation theory, and wave stability problems ( Bridget 1 1997al) . iBridgeJ 
1 1997lj) ). Bridges and ReicH ( 200d) develop multi-sy mplectic difference schemes. M ulti- 
symplectic Hamiltoni a n systems have bee n studied by Marsden and Shkoller ( 19991) and 
Bridges et al. (2005). Webb et al. ( 20071) : 11^666 et al. I ( 200^ 2014d) discuss tr a velini 


waves in multi-fluid plasmas using a multi-symplectic formulation. \Holm et al\ (1199 


give an overview of Hamiltonian systems, semi-direct product Lie algebras and Euler- 
P oincare equations. 

Cotter et al. I ( 200?!) developed a multi-symplectic, Euler-Poincare formulation of fluid 


mechanics. They showed that multi-symplectic ideal fluid mechanics type systems are 
related to Clebsch variable formulations in which the Lagrange multipliers play the role 
of canonically conjugate momenta to the constrained variables. Thus, the Clebsch variable 
formulation involves a momentum map. 
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The main aim of the present paper is to correct and ex tend the analysis of the multi- 

(2014c). Section 2 describes the MHD 


symplectic MHD equations derived hv lWebb et al 


equations and the first law of thermodynamic s. An overview of th e key equations of the 
multi-symplectic MHD equations obtained by Webb et al. ( 2014cll is given in Section 3. 
We correct some typographical errors and present a more consistent description of the 
Cartan-Poincare form for the system (section 4). The Cartan-Poincare form is used to 
obtain a set of differential forms representing the system using Cartan’s geometric theory 
of partial differential equations. 

Section 5 concludes with a summary and discussion. 


2. The Model 

The magnetohydrodynamic equations are: 


dp 


+ V-(/9u) = 0, 


^W + V. 


/9UU 


bb 


= 0 , 


dB 


+ U-VS = 0, 


— V X (u X B) + uV-B = 0, 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 


where p, u, p, S and B are the gas density, fluid velocity, pressure, specific entropy, and 
magnetic induction B respectively, and I is the unit 3x3 dyadic. The gas pressure p = 
p{p, S) is a function of the density p and entropy S, and p is the magnetic permeability. 
Equations are the mass and momentum conservation laws, (12.3|) is the entropy 

advection equation and (12.4p is Faraday’s equation in the MHD limit. In classical MHD, 
(ITTD-(ITD are supplemented by Gauss’ law: 


V-B = 0. (2.5) 

which implies the non-existence of magnetic monopoles. The MHD equations are closed 
by specifying an equation of state for the gas U = U{p, S) where U {p, S) is the internal 
energy per unit mass. The first law of thermodynamics has the form: 

TdS = dQ = dU + pdr where r = -, (2-6) 

P 

where t = 1/p is the specific volume. Using the internal energy per unit volume e = pU 
instead of U, (ITBl) may be written as: 

TdS = - {de — hdp) where h = ^ (2.7) 

P P 

is the enthalpy of the gas. Equation (IT71) gives the formulae: 

pT = es, h = Ep, p = pEp- e, (2.8) 

relating the temperature T, enthalpy h and pressure p to the internal energy density 
e(p, S). 
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3. Multi-Symplectic Appproach 

In this section we gi ve a brief overview of the multi-symplectic formulation of MHD 


hY \Wehh et al. (2014c). 


3.1. Clebsch Variable Action Principle 
Consider the MHD action (modified by constraints): 


J = 


d^x dtL, 


(3.1) 


where 


L = 


- e(p, S) 




-P 


+ r- 






| + V.(pu) 


u-V^ 


— - V X (u X B) + u(V-B( 


(3.2) 


In (13.21) the Lagrange multipliers </>, P, A and F ensure that the mass, entropy, Lin 
constraint and Faraday’s equation are satisfied. Note we do not set V-B = 0 in our 
analysis. It can be set equal to zero after all variational calculations for the MHD 
system are finished. However, it does modify the Lagrange multiplier equa t ion for F 

that ensures Faraday’s equation is s atisfied (see Morrison and Greend (Il980l). Mmrison _ 

19821) Morrison and Greene (Il982). Holm and Kuvershmidi i 1983a ) Holm and Kuvershmidi 


1983b ) and Chandre et al\ ( 2013h for further discussion of the MHD Poisson bracket). 


The Lagrangian in curly brackets equals the kinetic minus the potential energy (internal 
thermodynamic energy plus magnetic energy). The Lagrange multipliers </>, P, A, and F 
ensure that the mass, entropy, Lin constraint, Faraday equations are satisfied. 

Stationary point conditions for the action are 6J = 0. 6J/Su = 0 gives the Clebsch 
representation for M = pu: 


M = pu = pVp - PVS - AVp + B X (V X F) - FV-B. 


(3.3) 


Setting 6J/S(j), SJ/Sp, SJ/SX, SJ/ST equal to zero gives the mass, entropy advection, Lin 
constraint, and Faraday (magnetic flux conservation) constraint equations: 

pt + V-(pu) = 0, 

St + u-VS” = 0, 

Pt + u-Vp = 0, 

Bt - V X (u X B)+ u(V-B) = 0. (3.4) 

Similarly, setting 6J/dp, dJ/6S, SJ/6p, SJ/S'B equal to zero gives evolution equations 
for the Clebsch potentials p, /3, A and F as: 

'dp _ A 1 


dt 


+ u-\7p + -It — h = 0, 


(3.5) 


dB 

^ + V-(/3u)+pr = 0, 
f)\ 

^ + V.(Au)=0, 


f-„x(vxr) 


B 


■V(F-u) + — =0. 
Po 


(3.6) 

(3.7) 

(3.8) 





































4 G. M. Webb, J. F. McKenzie and G. P. Zank 

Equation ()3.5p is Bernoulli’s equation. The V(r-u) term in (j3.8ll is associated with V-B ^ 

0 . 


3.2. Multi-symplectic approach 


FTom \Webb et al. (j2014cr i. the MHD system was written in the multi-symplectic form: 


dz^ 

Q, 

ij 


dH 
dz^ ’ 


(3.9) 


where the variables z^ are defined as: 

z={u\u^u^p,S,pL,B\B^,B^T\T^,T\X,P,(l)f. (3.10) 

The variables: a:“ = {t, x, y, z) are the independent space and time variables, i.e. x^ = t, 
x^ = x, x^ = y, x^ = z. Below we use the notation = dz’^jdx'^. The constrained 
Lagrangian L in (13.211 is written in the form: 

L = Lo + E^“<- (3.11) 

S 

where 

-^0 = -£(p,5') - (3.12) 

is the unconstrained Lagrangian. 

The multi-symplectic Hamiltonian is given by the generalized Legendre transformation: 

B[z) = ^ -L = -Lo = - {^-pu^ - e{p, S) - ^ ■ (3.13) 


The fundamental one forms a;“ {a = 0,1, 2, 3) of the multi-symplectic system (13.91) are 
defined as: 


a;“ = L‘^{z)dzF 

Taking the exterior derivative of the a;“ gives the formulae 

= duj°‘ = lK,dz^ A dz\ 

2 ^ 

which, in turn, gives the formulae: 


K“ = 


dL° 


dz^ 


(3.14) 

(3.15) 


(3.16) 


for the skew symmetric matrices in (13.91) . 

For the variables z^ in (I3.10p the one forms a;“ are given by: 

uj^ =4>dp jddS Xdp F-dB, 

w* =K iPdS + Xdp F (j)dp) + pcj^dK + F-Bdu* - B^T-du) + u^r-dB), (3.17) 


where 1 ^ i ^ 3. 

The pullback conservation laws: 
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follow from the pullback of the identities: 

^ = d - Hiz)} = dL, 


= dL = ^dx^. (3.19) 

(see proposition 4.1 of liye66 et al. (2014c) for more detail). 

The pullback equation k°‘^ = 0 where k“ = do;" gives rise to the symplecticity or 
phase-space conservation laws (structural conservation laws): 




These conservation laws can also be written in the form: 


DpG-, - D^Gp = (K“ = 0, 


(3.20) 


(3.21) 


i.e., the symplecticity conservation laws (I3.2ip are compatibility conditions for the pull¬ 
back conservation laws (13.181) . 

The pullback conservation law (I3.18|) for /3 = 0 ( Wehh et al. I ( 2014c ). equation (5.39)) 
reduces to the energy conservation law: 

B^\ 


Go = - 




pu^ + e(p, S) + —j + V-[ pu 


1 


h 


E X B 


= 0, (3.22) 


where h = {e -\- p)/p is the enthalpy of the gas, E = —u x B is the electric field, and 
E X B//i is the Poynting flux. Similarly, the pullback conservation laws (I3.18|) for ,5 = i 
(1 ^ f ^ 3) give rise to the MHD momentum conservation equation: 


g- = -4w 


pu (g) U - 




B 


2p J 


= 0 , 


(3.23) 


Webb et al. lj2014c[ ). equation (5.41)). 


4. Extensions, comments and corrections 

The symplecticity conservation laws (I3.20I) - (I3.21I) have a generalized curl form. Con¬ 
sider the symplecticity laws (13.211) for 1 < i, fc ^ 3, namely: 

= AGfe - DkG, = 0. (4.1) 

Introduce the dual of the tensor defined as: 

EP = = -(V X Gf, (4.2) 

where V x G is the spatial curl of G. Taking into account the momentum conservation 
law (13.231) for G, (14.21) reduces to: 


V X G = - 


— V X M -f V X 
ot 


V- I M ® u — 


B (g) B 


= 0 , 


(4.3) 


where 

M = pu, (4.4) 

is the momentum density or mass flux M of the MHD fluid. Note there is no contribution 
from the magnetic pressure (B^/(2p)) and gas pressure (p) gradient force terms in (14.3|) 
because V x V(p -I- B'^/{2p)) = 0 when one takes the curl of the momentum equation 
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(13.231) . The evolution of V x M in (14.311 is thus determined by the inertia and magnetic 
tension components: 

(4.5) 


M( 


> u — 




of the MHD stress-energy tensor. This suggests that (14.31) describes Alfvenic type distur¬ 
bances, in which both fluid spin and magnetic tension forces are part of the dynamics. 
Equation ()4.3p can also be expressed in the conservation law form: 




V X M -I- Vs 


V X Mu^ - 




= 0 . 


(4.6) 


Pressure gradient forces play no role in the vorticity-like conservation laws (ig and 

(Irel) . 

The symplecticity conservation law (14.311 is different than that obtained by taking the 
curl of Euler momentum equation in the form du/dt = F where F is the net force on 
the flu id element, to obtain an equation for th e ev olution of the fluid vo rticity a> = 


V X u. Webb et al. (2014a) Webb et ai 1 2014bll and llTe66 and An^ ( 2015l l obtained a 


conservation law in ideal fluid mechanics (i.e. for B = 0) for the generalized vorticity 

rj = u; -I- Vr X VS where u; = V x u, (4-7) 

is the fluid vorticity and r satisfies the equation: 

dr f d 


dt 


dt 


+ u-V r = —T. 


(4.8) 


Here r = jd/p where jd is the Cleb yh potential that ensures dS/dt = 0 in the Eulerian, 
Clebsch variational approach (e.e:. \ Zakharov and Kuznetsoi (Il997 1) and d/dt = d/dt + 


u-V is t he L agrangian time de r ivativ e following the flow. I Webb et al\ (l2014alll lTe&& et al. 


( 2014bll and Webb and Ancd ( 2015ll show that for fluid dynamics (B = 0) the modified 
vorticity flux fl-dS is advected or Lie dragged with the flow, i.e.. 


-(a.ds) = 


dfl 

— -V X (u X H) -I- u(V-ri) 

at 


■dS = 0. 


(4.9) 


The conservation law 621) and the associated conservation law for the modihed fluid 
helicity u-fi are nonlocal conservation laws that depend on the nonlocal var iable r = 

— fg T( x., t)dt where the integ ration is with respect to the Lagrangian time t (e.g. l Webb et al. 


( 2014all : [tTe66 et ^ ( 2014bl ll. Conservation law (14.91) in fluid dynamics is analogous to 


Faraday’s equation in MHD. Note that V-il = 0 in (14.91) . 

There is another symplecticity conservation law obtained from (13.211) for the case /3 = 0 
and 7 = 1,2,3. In that case (I3.21|) reduces to: 


m 


G - VGo = 0, 


(4.10) 


where G = 0 is the momentum equation (13.231) and Gq = 0 is the energy conservation 
equation (13.221) . 

The general form of t he symplecti c ity eq uations for MHD using Eulerian Clebsch 


potentials were given in \ Webb et al. (|2014cll [equations (5.44) et seq. of that paper]. 
There were some typographical errors in t he flux Fjy, indicat ed below. The general 


symplecticity conservation laws obtained hv lWebb et al. \ ( 2014cll have the form: 




(4.11) 
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where 


and 


pO _ 

^ab — 


di^,p) , d{p,s) , d{x,p) , d{rs,B^) 


+ 


+ 


d{x°-,x^) d{x°-,x^) d{x°-,x^) 


+ 


d{x°‘,x’^) ’ 


(4.12) 


77'fc _ 

^ab — ~ 


d{pu^,(j)) d{l3u^,S) d{Xu^,p) 

d{x°-,x^) d[x°-,x^) d{x‘‘,x^) 

d{TsB‘^,u'^) 5(r,S'=,u") 


d{x°-,x^ 




d{x°-,x^ 


(4.13) 


In (I4.13|) 1 ^ fc ^ 3. The fou r th term on the right handside of (I4.13P was missed in 
equation (5.48) in Webb et al. I (l2014cli . Also in (14.131) we used the identity: 


d{(jm^,p) d{p(j),u^) d{pu^,(j)) 


(4.14) 


9(x“,x^) 9(x“,x^) d{x°‘,x^)’ 

to simplify (5.48) of Webb et al. I ( ^14cll . The derivation of the symplecticity conservation 
laws (14.31) and (14.101) using the general symplecticity laws (14.111) and using (I3.3I) - (I3.8I) is 
a non-trivial algebraic exercise. 

4.1. Differential forms approach 

Proposition 4.3 in Webb et al. 1 1 2014cl l contains flaws. A consistent approach to the multi- 
sy mplect i c equa t ions using differential f orms for ID Lagrangian gas dynamics was given 
by I WebA ( 2015h . Webb and Ancd ( 2015 ) have given the corresponding theory for multi¬ 
dimensional, ideal, compressible, Lagrangian gas dynamics. Below we us e differential 


forms t o describe the Eulerian, Clebsch variable MHD variational principle of lITe66 et al. 

touch . 

Proposition 4.1. The multi-symplectic system hS.Wl is a stationary point of the action: 

J = j ^*(0) = J LdV, (4.15) 

where ^p*{Q) is the pullback of the differential form 0 given below, namely: 


0 =a;“ A - iLdP, uj°^ = L'^dzf 
dV =dt f\dx Ady f\ dz, dx°‘ = da-idV = (—l)“dxo A ... dx°‘~^ A dx 




A dx", 
(4.16) 

where we use the notation (x°,x^,x^,x^) = {t,x,y,z), and L is the constrained La¬ 
grangian 

Proof.The pullback of the form 0 is given by: 

^*(0) =fj* (L^dz^ A dxa - HdV) 


=Lf^dx^ A dx°‘ - HdV. 

^ dx^ 

However, 

dx® A dxa = dx® A (—l)“dx° ... A dx““^ A dx°‘~^^ ... A dx" = (—l)^“d® dE. 
Thus, 

dz^ 


riQ)={L-^-H]dV^LdV, 


(4.17) 

(4.18) 

(4.19) 


where L is the multi-symplectic Lagrangian (13.101) . 
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The stationary point conditions, 5JjSz'' = 0, give the Euler-Lagrange equations: 




dL 


d / dL 


dxi \ dz'-. 


= K“^-^ = n 
dx^^ dz^ ’ 


which is the multi-symplectic system (|3.9p (see also lHvd^ (2QQ5)). 
Proposition 4.2. Gonsider the variational functional: 


(4.20) 

□ 


G[i 2 ] = / n, 

J M 


(4.21) 

(4.22) 


where 

Q = dQ = duj°‘ A dxa — dH A dV. 

and M is a region of the jet space (fiber bundle space) with boundary dM, in which the z®* 
are taken as independent of the base variables x°‘ (a = 0,1, 2, 3). Gonsider the variational 
principle: 


where 


dCP = [ Cv (LI) = 0, 

JM 


r - ^ - v^ ^ 


(4.23) 


(4.24) 


is the Lie derivative with respect to the arbitrary, but smooth vector field V. The varia¬ 
tional equation (5G[r2] = 0 reduces to: 


(5G[12] = / = 0, 


IdM 


where the forms ftp are given by the formulae: 

d 


dH 




(4.25) 


(4.26) 


(\ ^ p ^ N). Because the are arbitrary smooth functions of the z®, the variational 
principle <5G[f2] = 0 implies: 

Pp = 0, l^pi^N. (4.27) 

The pullback of the forms {dp} to the base manifold gives the equations: 

dz^ dH' 


/?p= K 


dx^ dzP 


dV = 0. 


(4.28) 


Thus, the sectioned forms ftp vanish on the solution manifold of the multi-symplectic 
system iTOI) . and the {dp} can be used as a basis of Gartan forms describing the system 

Proof.The proof is essentially the same as that given bv llTe6tl ( 2015 1 for the case of ID 
gas dynamics. A critical component of the proof is the use of Cartan’s magic formula: 


£v (fl) = VjrfD + d (Vjfl) = d (Vjfl), (4.29) 

where we used the facts Ll = d0 and dLl = dd0 = 0. Using (14.291) and Stokes theorem, 
(|4.23p reduces to: 

dG[f2] = [ dCVM) = f VM = [ VP( ) = [ VP dp = 0, (4.30) 

^ ' Jm ^ ’ JdM JdM \d^P JJom '^ ^ ’ 
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which verifies (j4.25|) . The formula (14.2611 for (3p is obtained by using (j3.15|) for c?a;“ and 
(j4.22l) for n, to obtain: 




dzP 


dzP' 


1 dH 
-K“ (iz* A dz^ A dxa — -z — dz^" A dV 

2 dz^ 


(4.31) 


Using the skew symmetry of K“ and dz* A dz ^, (I4.31|) reduces to the expression (I4.26|) 
for /3p. This completes the proof. □ 

4.2. The differential forms ftp 

The differential forms ftp = dzpzil in (|4.26|) may be used to represent the MHD system 
described by the Clebsch variable variational principle. The d ependent variables z are 
listed in (|3.10|) . In the time evolution variational principle (e.s. \ Zakharov and Kuznetsoi 
(|l997l) . the fluid velocity u is expressed in terms of the Clebsch potentials, and is elimi¬ 
nated from the Hamiltonian density id = {\/2)pv?£{p, S)B'^/{2yL) and {p,4>), iS,ft), 
(/i,A), (B,r) are canonically conjugate pairs in the canonical Poisson bracket. We use 
the notation: 

(3^'=dz.M, (4.32) 

where the Cartan Poincare form in (j4.22l) has the form: 


OH 

fl = A dxo -b A dxk — t: — dz^ A dU, 

ozP 


(4.33) 


where the differential forms and uj^ are listed in (13.171) . From (14.321) and (14.331) 

dH 


= dziz{ duj° A dxo + dw^ A dxk — Tr~dz^ A dU 

ozP 


(4.34) 


Using (13.171) and (14.341) we obtain: 

/3“ = {ftdS -b Xdp — pdf) — B’^dTs) A dxi -b {VidB^ -b B^dTi) A dxk + pu^dV, (4.35) 

for the differential forms associated with u. Using the identity 

dx°' A dxi = SfdV, (4.36) 

the sectioned form equation (3^ = 0 yields the expression: 

pu = pVf) - fVS - AV/r -b B-(Vr)^ - B-VP - PV-B, (4.37) 

which is equivalent to the Clebsch expansion for the mass flux pu given in (13.31) . 

The differential form ft^ is given by: 


OH 

pP = = (9pjda;°) A dxo + {dpzduj^) A dx^ - 


dp 


Using (13.141) for H, we obtain: 

dH 


dp 


= - \ - £, 




where h= {e -\-p)/p is the gas enthalpy. Substituting (14.391) in (14.381) gives 
pP = — {df) A dip + u^df) A dxk) + ^ 


PP=- 


df) 

dt 


- -u^ -h 


dU, 


(4.38) 

(4.39) 

(4.40) 

(4.41) 
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for the differential form pP and for the sectioned form pP. Note that /3^ = 0 is equivalent 
to Bernoulli’s equation (|3.5p . Similarly, we obtain: 

=ds-ifl = —dp A dxo — d{Pu^) A dxk — pTdV, 


(^+V-(/3u) + pt) dV. 

The equation P^ = 0 corresponds to (13.61) for p. 

Following the above procedure we obtain the equations: 

pP =df^zfl = — [dA A dxQ + d{Xu^) A dxk\ , 


(4.42) 

(4.43) 


B'- 


dVi A dxo + u dVi A dxk + F^du® A dxi H- dV 




P^* =dB^ A dxo + {B^du^ + u’^dB^ - B^dK) A dxk, 
P^ =9a-iA1 = {dp, A dxo + u^dp A dxk) , 

P^ =dpzGl = dS A dxo + u^dS A dxk, 

P'^ =drj)zVl = dp A dxo + {u^dp + pdu^) A dxk- 


(4.44) 


The pullback of the above equations, i.e. P^” = 0, gives the evolution equations for 
{\,F,'B, p, S, p) listed in (I3.4p - (I3.8I) . The differential form equations P^^ = 0 thus repre¬ 
sent the partial differential equation system (1^ . 

It is not obvious that the system of forms {P^’’} above is a closed ideal. A check on 
the closure of the forms for the case of non-barotropic, ID gas dynamics (i.e. B = 0) 
indicates that the ideal of forms I = {P"^, pP, P^, P^, P‘^} can be closed by adjoining the 
form d/3“. The ideal I is closed for the case of a barotropic gas. The Cartan approach to 
Lie symmetries requires that I is a closed ideal (e.g. Harrison and Estabrook (1971)). The 
ideal is closed if dPi = Ca A P ^^, w here the are forms. It should be noted that the ideal 


of forms obtained by WebA ( 2015 1 for ID, Lagrangian, multi-symplectic gas dynamics is 
closed. The ideal of forms for multi-dimensional, Lagrangian. compressible gas dynamics 
obtained from the Cartan-Poincare form is also a closed ideal (Webb and Anco (2015)). 
This suggests that the ideal of forms using the Clebsch variable description has a more 
complicated structure than the set of forms that arise in the Lagrangian variational 
approach. In general, this question requires more work and lies beyond the scope of the 
present paper. 


5. Summary and Concluding Remarks 

In this paper, we have corrected some flaws in the paper bv I Webb et a,L I ( 2014c 1 on 
multi-symplectic MHD. This includes the correction of some typographic errors in the 
expressions for the conserved fluxes in the symplecticity conservation laws written in 
terms of the Clebsch potentials in (I4.11I) - (I4.13I) and in particular expression (14.131) for 
has been corrected. We pointed out that the vorticity-type symplecticity conserva¬ 
tion law is equivalent to taking the curl of the MHD momentum equation in the form 
(gH). This conservation law has conserved density V x M where M = pu is the momen¬ 
tum density of the MHD fluid. Pressure gradient forces due to the combined gas and 
magnetic pressures (i.e. V(p + ps), where pb = B^/(2p) is the magnetic pressure and 
p is the gas pressure), play no role in this conservation law. The conservation law (14. 6|) 
describes fluid motions with spin and rotation (e.g; as in an Alfven wave) . This symplec¬ 
ticity conservation law is different than that obtained by Webb and Anec ( 2015h in ideal 
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Lagrangian fluid mechanics. In the latter paper the vorticity symplecticity conservation 
law involved O = uJ + VrxViS” where u; = V x u is the fluid vorticity and r = /3/p where 
P is the Clebsch potential that enforces entropy conservation following the flow. This is 
a nonlocal conservation law as it involves the nonlocal variable r. It has the same form 
as Faraday’s equation in MHD except that B is replaced by $7. A further symplecticity 
conservation law is due to a compatibility condition which requires that time derivative 
of the the momentum conservation equation minus the gradient of the energy conserva¬ 
tion equation is zero. We provided a more consistent treatme nt of variational prin ciples 


associated with the Cartan-Poincare form than that given in I Webb et al. ( 2014(J) . and 


its relation to Cartan’s geome tric theory of partial d i fferen tial equations using differential 
forms (Sections 4.1 and A.2). \Cendra a,nd Ca,vriot\ (2013) suggest that the formulation 
of multi-symplectic systems using exterior differential forms, is a natural approach in 
taking into account integrability conditions in the theory. 
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